Light funneling mechanism explained by magneto-electric interference 
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We investigate the mechanisms involved in the funneling of the optical energy into sub- wavelength 
grooves etched on a metallic surface. The key phenomenon is unveiled thanks to the decomposition 
of the electromagnetic field into its propagative and evanescent parts. We unambiguously show that 
the funneling is not due to plasmonic waves flowing toward the grooves, but rather to the magneto- 
electric interference of the incident wave with the evanescent field, this field being mainly due to 
the resonant wave escaping from the groove. 



Plasmonics, as the science of the efficient coupling of 
photons with free electron gas oscillation modes at the 
surface of metals, appears as an inescapable solution for 
the design and realization of optical nano antennas pQ. 
Numerous cutting edge applications are based on nano- 
antennas like biosensing [2], gas sensing [3], photovoltaic 
[4] or infrared photodetection [5] which exploit the in- 
tense local electromagnetic field in a confined volume [Bi- 
ll]. Now, the specific matter of total photon harvesting 
at the nanometric scale, i.e. designing an antenna able 
to couple all the incident optical power with a nanoab- 
sorber, remains challenging [TJ [THS]- The natural two- 
step antenna sequence (collection of light, then concen- 
tration) has been extensively studied in structures made 
of a metallic subwavelength grating surrounding a tar- 
get [6] [T0UT4] . The underlying mechanism involves SPP 
excitation (collection) and propagation (concentration) 
along the grating until the coupling with the target. Such 
structures, though, are designed to collect light at a spe- 
cific incidence angle, which is obviously a strong practical 
limitation. 

In contrast, quasi-isotropic perfect transmission is ob- 
tained through very narrow slits drilled in a metallic 
membrane [TSJ [16] . This perfect transmission is success- 
fully explained by a localized Fabry-Perot resonance in 
the slits [17] . However the funneling, namely the mecha- 
nism responsible for the redirection and subsequent con- 
centration of the whole incident energy flow, from the 
surface toward the tiny aperture of the slits, remains un- 
clear. Yet, a pictorial model of the underlying physics is 
of a key importance for the design of efficient nanoanten- 
nas. 

Such a model is given by the energetic point of view 
[19]: Poynting vector streamlines distinctly show that 
the incident flow bends when reaching the metal sur- 
face, and then propagates along the interface toward the 
slits. This fits with the intuitive explanation, inspired 
by the SPP excitation process, that plasmonic waves 
drive the funneling sequence [20 . Furthermore, quasi- 
cylindric waves were recently identified as the dominant 
short-range propagation process of the field amplitude 
along the surface of the grating J2TJ [22]. Nevertheless, 



even if the evanescent waves are naturally assumed to 
concentrate the energy toward the apertures of the slits, 
no specific study of the light funneling has so far been 
carried out to our knowledge. 

In this letter, we definitely unveil the funneling pro- 
cess, and highlight the unexpectedly limited role of the 
evanescent waves alone. Firstly, by analysing the par- 
ticular case of a groove (for which experimental study 
confirm theoretical predictions [H]), we show that the 
evanescent waves do not carry any energy through the 
apertures: they simply redistribute it over the metal sur- 
face. Instead, we identify the magneto- electric interfer- 
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FIG. 1. (a) Poynting vector streamlines are drawn on a 
grating of grooves of height h = 640 nm, width w = 56 nm 
and period d — 2 /im. The incident energy is funneled inside 
the grooves, (b) Vector streamlines inside a groove. The 
dissipation of the energy is computed in the metallic region, it 
clearly appears that this dissipation occurs on the sidewalls of 
the grooves (orange volume), see [18 for detailed field maps, 
(c) Reflectivity spectra for various values of the grating period 
d — 1.618 /im, d = 2 /im and d — 2.472 /im. The geometry of 
the groove is the same as before. The period has almost no 
influence on the resonance since it is due to the Fabry-Perot 
resonator inside the grooves. 
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ence (MEI) [23] of the incident wave with the evanescent 
field as the main mechanism of the funneling sequence. 
We then use a single interface analysis [24 to generalize 
our result to sub- wavelength apertures with no resonator 
behind. MEI also explains the broadband extraordinary 
transmission due to plasmonic Brewster angle that was 
recently published by Alu [25] (see supplemental material 

Let us consider an infrared light at A f = 4 /im incident 
onto nanometric sized grooves periodically drilled into a 
gold surface. Figure[l](a) shows the geometry of the grat- 
ing (width w = 56 nm, height h = 640 nm and period d) . 
The period is chosen so that there is no diffracted wave 
for all angles of incidence (hence d < A//2). The light is 
TM-polarized (transverse magnetic) and incident with an 
angle 0. The dielectric function of gold is computed from 
the Drude model s(X) = 1 — [(X p /X + ij) A p /A] _1 which 
is suited to the infrared spectral range for X p = 161 nm 
and 7 = 0.0077 [26 . The electromagnetic analysis of 
this structure is done using a B-spline method [27] , which 
can perform fast and exact computation of Maxwell equa- 
tions. The Poynting- vector streamlines show how the en- 
ergy flow is funneled toward the apertures, and dissipated 
mainly on the sidewalls of the grooves. The reflectivity 
of the grating is plotted in Fig. [I] (c) at normal inci- 
dence for d = Xf /2 and for a random value d = 1.618/im. 
Although the grating is structured on a tiny portion of 
its surface (less than 3%), it exhibits a resonance with a 
total absorption at normal incidence. We should high- 
light that Xf depends only slightly on the period d. The 
absorption remains nearly total even for large incidence 
angles (0 < 50°), which is adapted to light collecting 
systems. 

In order to address the funneling mechanism, we con- 
sider the electromagnetic field in the air, and we split it 
into three terms. The magnetic field is expressed: 

Htotai = Hi + H r + H ei (1) 

where Hi is the incoming wave, H r is the reflected wave, 
and H e is the sum of the diffracted evanescent waves. 
Similar definitions can be given for the electric field com- 
ponents. In the rest of this letter, E x H stands for 
the mean time average value of the vectorial product 
and is practically computed from complex amplitudes as 
^iR(E x H*). Thanks to the decomposition of Eq. |l]), 
the Poynting vector can be expressed as the sum of six 
terms: 

S — Si -\- S e i -\- S r -\- S er -\- S e -h Si r) 

(2) 

with Si = Ei x Hi, S ei = E e x Hi + Ei x H e , S r = 
E r x H r , S er — E e x H r -\- E r x H ei S e = E e x H e , and 
Si r — Ei x H r + E r xHi. The terms Si and S r are respec- 
tively the incident and the reflected fluxes of the plane 
wave. The term S e corresponds to the energy carried by 
the evanescent waves. The term S e i corresponds to the 



MEI between the evanescent and the incident fields. It is 
easy to prove that the six terms of equation Q are flux 
conservative (null divergence) thus each of them can be 
considered to be an independent energy flux vector in the 
air. In order to simplify the discussion, we first consider 
an optimized device at the resonance wavelength. So we 
have no reflected wave, the fields H r and E r are null and 
the Poynting vector can be expressed as S = Si+S e i + S e . 

The Poynting- vector streamlines for Si, S e i and S are 
plotted at two angles of incidence in Fig. Q so that the 
flux of energy between two lines is constant. 

At normal incidence, as expected for a propagative 



4,0° V 30 ° 




0.5 1 1.5 0.5 1 1.5 

x(|jm) x(|jm) 



FIG. 2. Poynting- vector streamlines in one period of the slit 
grating for two angles of incidence — ° (left column) and 
— 30 ° (right column) at A = 4000 nm. Streamlines of the 
incident wave are shown on (a) and (b). Streamlines of the in- 
terference between the incident wave and the evanescent field 
are shown on (c) and (d). The energy flux of the evanescent 
waves is negligible in this structure for ^ = 0°; refer to Fig. 
[3] (a) for an illustration at — 30 . Streamlines of the total 
Poynting vector are shown on (e) and (f). In both cases, the 
incident energy is funneled inside the groove where it is fully 
absorbed inside the metal. 
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plane wave in air, the lines for incident flux Si are equidis- 
tant and in the propagation direction. The MEI S e i lines 
are coming from the surface and are converging on the 
groove. On the metallic surface, they compensate for the 
flux of the incident plane wave and funnel it inside the 
groove. By drawing lines perpendicular to the Poynting 
streamlines on Fig. [2](c), and taking in account that the 
predominant term in S e i is E e x Hi, one can deduce that 
the evanescent wave shape is quasi-cylindrical [21j [22] . 
The evanescent flux S e (not shown) carries energy 1000 
times weaker and does not play an active role in the fun- 
neling for this structure at normal incidence. In Fig. [2] 
(e), the total flux of energy S is shown to funnel into the 
groove in the near- field region [z < 500 nm). Eventu- 
ally all the incident flux is dissipated, mostly inside the 
groove. 

For an incidence of 30 ° , the MEI S e i is still funneling 
the energy towards the slits (Fig. [5](d)). But it no longer 
compensates for the incident flux which lines are equidis- 
tant. Indeed, in Fig. [2] (d), there are more lines going 
out from the metal surface on the left (10 lines) than on 
the right (6 lines). Nonetheless, Fig. |2](f) shows that the 
incident energy gets funneled into the groove despite this 
assymetry. In fact, at oblique incidence, the evanescent 
field carries an energy flux S e which is no longer negli- 
gible, as shown in Fig. [3j the energy is redirected from 
the right side of the groove to the left side. This redistri- 
bution of the energy compensates for the dissymmetry of 
S e i which appears in Fig. [5](d). In conclusion S e plays no 
role in the funneling, but helps to redistribute energy over 
the grating. This incidentally invalidates the hypotheti- 
cal role of plasmonic waves, which are evanescent waves, 
in the funneling mechanism. As a general comment, it is 
interesting to point out that the MEI process is known to 
be responsible of the optical tunnelling (frustrated total 
reflection) [28 . Our study unveils its key role in a larger 
spectrum of near-field energy transfer phenomena. 
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FIG. 3. Poynting- vector streamlines of the evanescent field 
for two angles of incidence (a) 30 ° and (b) 50 ° at A = 
4000 nm. S e does not play a role of funneling but redistributes 
the energy in the grating. 

We now aim to describe the origin of the evanescent 



field involved in the funneling process. The subwave- 
length groove behaves as a Fabry-Perot resonator. As 
in ref. [24], we consider the isolated single interface in 
two configurations: one where the incident field is a unit- 
amplitude plane wave in air; the other where the incident 
field is a wave coming from the bottom of the groove, 
and is unit-amplitude at the interface. In the first case 
(Fig. [4] (a)) the reflected wave has an amplitude p\ and 
an evanescent field 771. Due to the low aperture ratio 
w/d, we get \pi\ < 1 and 1 771 1 <C 1. In the second case 
(Fig. [4] (b)) the transmitted wave (which corresponds 
to the reflected wave above) has an amplitude T2 and 
the evanescent field in the air is written 772. If a unit- 
amplitude wave is defined as taking the value H y = 1 at 
the center of a groove entrance, then computation shows 
that evanescent amplitudes 771 and 772 are nearly equal. 
Due to the low aperture ratio w/d, we get |t2 | <C 1. 
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FIG. 4. Isolated single interface analysis of the metallic grat- 
ing optimized to absorb all the incident light, (a) Unit plane 
wave from air: pi is the amplitude of the reflected plane wave, 
771 is the vector of evanescent field amplitudes, (b) Unit modal 
wave from bottom of grooves: T2 is the amplitude of the plane 
wave escaping in the air, 772 is the vector of evanescent field 
amplitudes, (c) Same as (b), but with the modal wave having 
the amplitude A. (d) Superposition of (a) and (c), showing 
the field amplitudes in the grating excitated by an unit plane 
wave: the reflectivity is null and the evanescent field is dom- 
inated by the term escaping from the resonator. 

At the resonance, the wave coming from the bottom 
of the groove has an amplitude A — —r^ 1 pi so that all 
the amplitudes at sake in Fig. [4] (b) are multiplied by 
the factor A, which leads to Fig. [4] (c). The response 
of the grating excited by a unit-amplitude plane wave 
at this resonance is thus given by the superposition of 
amplitudes of Fig. [4] (a) and (c), shown on Fig. [4] (d). 
This leads to the expected null-amplitude reflected wave. 
Moreover, the resulting evanescent field is expressed as 
771 + 772 A ~ 772 A because \rji\ ~ |t7 2 | and A > 1 (e.g., for 
the grating described in Fig. [I] one computes \A\ ~ 11). 
To summarize, the wave built inside the grooves escapes 
in the air as both a propagating plane wave r^A and an 
evanescent field 772 A (see Fig. [4] (c)). The propagat- 
ing plane wave interferes destructively with the directly 
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reflected wave (J2A + pi = 0) leading to the null reflec- 
tion, and the evanescent field interferes with the incom- 
ing plane wave to funnel the energy into the groove. The 
two effects are of course not independent: we have no 
reflection because the energy is nearly fully collected in 
the groove. As a generalization path, it is interesting to 
note that this analysis still holds whatever the structure 
etched behind the aperture on the metal surface (rectan- 
gular slits, or more complex shapes such as in Ref.[29 ). 

Finally, we want to highlight that the key role played 
by the MEI of propagative waves with the evanescent 
field in the funneling mechanism is not limited to reso- 
nant structures. Let us thus consider the non-resonant 
case of a grating made of infinitely deep grooves, i.e. the 
simple isolated interface illustrated in Fig. |4ja), with a 
reflected wave of amplitude p\. For an incidence angle 
of 30°, we get |pi| 2 = 0.83: about 17% of the incident 
energy enters the grooves. This value is much larger than 
the aperture ratio w/d = 2.8%, thus there appears to be 
funneling. Now, if we compare this to the situation of 
Fig. [4^d), the evanescent field of Fig. [4ja) is lowered by 
a factor \rji -\- rj2A\/\rji\ ~ \A\ ~ 11. We therefore com- 
pute that the interference of the incident wave with the 
evanescent field gives a funneling of about 1/\A\ ~ 9%. 
The missing 8% stems from the interference of the evanes- 
cent field with the reflected wave pi , which paradoxically 
contributes to the funneling of energy inside the grooves. 
The key to the paradox lies here: firstly, the main con- 
tribution to the interference is E e x H r (E r x H e is much 
smaller, at least at the interface level) and, secondly, the 
magnetic field H r has the sign of Hi due to the metallic 
reflection. 

In conclusion, we have unveiled the funneling mecha- 
nism of incident light in very narrow grooves etched on 
a metal surface. It originates from the magneto-electric 
interference between the incident wave and the evanes- 
cent field, in both resonant and non resonant situations. 
Furthermore, this result has been generalized to any sub- 
wavelength aperture etched on a metal surface (whatever 
the structure behind it) thanks to a single interface anal- 
ysis. In the resonant case, the evanescent field escaping 
from the apertures can lead to the full harvesting of in- 
cident photons for a broad range of incidences. From 
a practical point of view, this approach opens a new av- 
enue for the design of electromagnetic resonant antennas, 
based on the tailoring of the escaping evanescent field. 

Eventually, we have shown that evanescent waves prop- 
agating along the interface do not carry any energy 
through the apertures. This clearly demonstrates that 
the funneling is not mediated by plasmon waves at the 
surface. 

This work was partially supported by the ANTARES 
Carnot project. 
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SUPPLEMENTAL MATERIAL 

Supplementary information for the article Light fun- 
neling mechanism explained by magneto- electric interfer- 
ence (MEI). It includes the following items: 

• Broadband extraordinary transmission analysis by 
MEL 

• Experimental results on groove gratings described 
in the article. 



• Field maps. 

• Scaling properties. 
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Broadband extraordinary transmission analysis by 
MEI 



The total funneling in the case of broadband extraor- 
dinary transmission described by Alu et al. [25] is illus- 



trated on figure |S2| and fully explained by MEI as shown 
on figure [ST] In constrast with fig. 3, the evanescent 
fields (amongst these evanescent fields are the surface 
plasmon contribution), carry here strictly no energy. As 
a result, in both normal and Brewster incidences, MEI 
alone redirects all the incident energy towards the slits. 



FIG. S2. Poynting vector streamlines in the broadband ex- 
traordinary transmission with parameters of Ref. [25] (wave- 
length A = 759 nm, angle of incidence = 68.60°, metal 
permittivity e m = —24.7+ 1.447, slit width w = 24 nm, pe- 
riod d — 96 nm, height h = 200 nm). Incident wave arrive 
from the top left. There is strictly no resonance inside the 
slit and no reflexion at all. The energy flux between two lines 
is constant, equal to 1/8 of the incident energy flux on one 
period. Lines need to be cut around a metal bloc (here at 
bottom of the right slit) because of the losses. 



Experimental results 




FIG. SI. Decomposition of incident energy flux of figure |S2] 
in two terms: the incident energy flux (straight green lines), 
and the magneto-electric interference (MEI, curved red lines). 
All other terms are null. Hence, the total funneling is fully 
explained by MEI energy collection from the surface to the 
slits. 



The theoretical computations presented in the article 
are confirmed by experimental results. To fabricate such 
a grating of high aspect ratio metallic grooves, we have 
developped a mold cast technique based on a gold electro- 
plating. The mold in GaAs is fabricated trhough induc- 
tively coupled plasma reactive ion etching (ICP-RIE). A 
scanning electron microscope (SEM) image of the GaAs 
walls is shown in figure [S3fc , they are 2/im high and 
150nm wide in a period d = 2.5/xm. A gold electroplating 
is then done on the GaAs mold which is eventually chemi- 
cally removed, leading to the grating of gold grooves (see 
SEM image in figure |S3}p). Measured reflectivity spec- 
tra of a TM incident wave at various angles (0 = 5°, 
6 = 30° and = 40°) are shown on figure Iggk They 



were obtained using a Fourier transform infrared spec- 
trometer and an home-made achromatic optical system 
[30] . As predicted, there is an almost total extinction 
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of light at A = 10.05/mi which is nearly independent of 
the incidence angle. The measurements are in excellent 
agreement with the computed spectra. The correspond- 
ing funneling behavior of the poynting energy flux are 



shown for these three incidence angles in figure S3 i-f. 
The reflected energy is neglected in these three cases. 





10 12 
Wavelength (|jm) 

(e) V 30° (f) 



^ 40° 




FIG. S3, (a) SEM image of the grating of GaAs walls which 
is used as a mold, (b) SEM image of the grating of gold slits 
(w — 150nm, h — 2/im and d — 2.5/im) obtained after an 
electroplating on the GaAs mold, (c) Reflectivity curves of 
the sample for various incident angles = 5°, = 30° and 
— 40°. The light is TM-polarized. The resonance peak 
wavelength is nearly independent of the incidence angle and 
leads to an almost total absorption. Poynting vector stream- 
lines at the resonance peak A = 10.05/im for an incidence of 
(d) 6 = 5°, (e) A = 30° and (f) 6 = 40°. The reflected energy 
is neglected. 



Field maps 

The magnetic and electric fields maps are shown in 
figure S4 a-b for the groove of width w = 56nm, depth 
h = 640nm and period d = 2/im at the resonance peak 



wavelength X p = 4/im. The incident wave is at normal 
incidence and in TM polarization. The magnetic field 
H y is concentrated at the bottom of the groove (z = 
— 640nm) while the electric field is concentrated at the 
aperture of the groove (z = Onm). 
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FIG. S4. Magnetic and electric fields ((a): \H y \ 2 , (b): \E\ 2 ) 
in the groove of width w = 56nm, depth h = 640nra and 
period d = 2/im at the resonance peak X p = 4/mi at normal 
incidence. 



Scaling properties 

The dependence of the resonance peak as a function 
of the groove's depth h is studied in the reflectivity map 
plotted on figure |S5| The width of the groove w = 56nm 
and the period of the grating d = 2/im are kept constant. 
The incident wave is TM-polarized and has an incidence 
of = 0°. There are three lines with total absorption 
of light appears. They correspond to the first three cav- 
ity modes inside the grooves. These resonances position 
as a function of the groove's depth is given by the law 
A m = 4/in e ff / (2m — 1) + <\> where m is an integer standing 
for the resonance order, n e ff is the effective index of the 
guided mode inside the groove and <ft is the wavelength 
shift induced by the reflection of the guided mode on the 
grating surface interface. At the first order, the effective 
index takes a very simple form: 



^eff 



1 + 5/% 



(3) 



where S is the skin depth of the metal [3T] . 

The dependence of the resonance peak as a function 
of the groove's width is studied in the reflectivity map 
plotted on figure S6 The depth of the groove h = 640nm 
and the period of the grating d = 2/mi are kept constant. 
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FIG. S5. Reflectivity map as a function of the groove's depth. 
Other parameters are kept constant (width w = 56nm and 
period d = 2/xm). The incident wave is TM-polarized and 
coming at normal incidence. 
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FIG. S6. Reflectivity map as a function of the groove's width. 
Other parameters are kept constant (depth h = 640nm and 
period d = 2/xm). The incident wave is TM-polarized and 
coming at normal incidence. 



The inverse dependance of the resonance peak with the 
groove's width is coherent with the previous law given in 
equation [3j The dependence of the resonance peak as a 
function of the grating's period is studied in the reflectiv- 
ity map plotted on figure S7 The depth h = 640nm and 
width w = 56nm of the groove are kept constant, which 
set the resonance peak at X p ~ 4/im for subwavelength 
period. The reflectivity map confirms this weak influence 
of the period value on the resonance peak. 
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FIG. S7. Reflectivity map as a function of the grating's 
period. The groove's width and depth are kept constant (w = 
56nm and h — 640nm). The incident wave is TM-polarized 
and coming at normal incidence. 



